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ABSTRACT 

We present several exact solutions of the eigenfrequency problem for torsional 
shear vibrations in homogeneous and non-homogeneous models of the neutron 
star crust governed by canonical equation of solid mechanics with a restoring 
force of Hookean elasticity. Particular attention is given to regime of large 
lengthscale nodeless axisymmetric differentially rotational oscillations which 
are treated in spherical polar coordinates reflecting real geometry of the neu- 
tron star crust. Highlighted is the distinction between analytic forms and 
numerical estimates of the frequency, computed as a function of multipole de- 
gree of nodeless torsional oscillations and fractional depth of the crust, caused 
by different boundary conditions imposed on the toroidal field of material dis- 
placements. The relevance of considered models to quasiperiodic oscillations, 
recently detected during the flare of SGR 1806-20 and SGR 1900-1-14, is dis- 
cussed. 
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1 INTRODUCTION 

Ever since the identification of pulsars witli neutron stars, tlie non-radial torsional sliear 
oscillations restored by bulk forces of different in physical nature internal stresses have 
been and still are among the most important issues in the study of the interconnection 
between the electromagnetic activity and asteroseismology of pulsars (e.g. Ruderman 1969, 
van Horn 1980; Hansen and Cioffi 1980; Schumaker and Torne 1983; McDermott, van Horn 
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& Hansen 1988, Strohmayer 1991; Duncan 1998; Bastrukov, Weber, & Podgainy, 1999; 
Yoshida, Lee 2002, Bastrukov, Podgainy, Yang & Weber, 2002; Bastrukov, Chang, Mi§icu 
§, Molodtsova & Podgainyi 2007). Recently, increasing interest in this domain of research 
has been prompted by the discovery of quasiperiodic oscillations (QPOs) on the lightcurve 
tail of SGR 1806-20 (Izrael et al 2005) and SGR 1900+14 (Strohmayer & Watts 2006). It 
has been suggested in these latter works that detected variability can be explained as set by 
quake induced torsional shear oscillations of neutron star. Different aspects of this discovery 
arc currently under intense theoretical investigation (Piro 2005, Glampedakis, Samuelsson & 
Andersson 2006, Watts & Reddy 2007, Lee 2007; Levin 2007; Sotani, Kokkotas, & Stergioulas 
2007; Samuelssen & Andersson 2007; Vavouhdis, Stavridis, Kokkotas, & Beyer, H. 2007). 
An extensive survey of observational data and theoretical works devoted to this issue can 
be found in (Israel 2007; Watts and Strohmayer 2007). 

Motivated by the above interest, we focus here on the mathematical physics of the eigen- 
frequency problem for the torsional shear oscillations governed by equation of Newtonian 
solid mechanics. The restoring force is the bulk force of Hookean elastic shear stresses. Em- 
phasis is laid on the boundary conditions which must be imposed on the toroidal field of 
material displacement at the edges of the seismogenic layer, that is, on the core-crust bound- 
ary and the star surface. Working from the homogeneous crust model we show that these 
boundary conditions substantially affect the asymptotic spectral formulae for the frequency 
of torsional shear oscillations. 

The plan of this paper is as follows. In Sec. 2, a brief outline is given of the governing 
elastodynamical equation for a standard core-crust model of a quaking neutron star with 
homogeneous crust. The general solution of the Helmholtz equation for toroidal field of 
material displacements describing the standing-wave regime of torsional oscillations in the 
spherical polar coordinates related to the real geometry of the neutron star crust is presented. 
Also, we obtain here the frequency spectrum for the nodeless global torsional elastic mode in 
the entire volume of the fiducial homogeneous solid star model which was also obtained in our 
previous investigation but by use of Rayleigh's energy method. In Sec.3, we derive two exact 
dispersion equations for the standing- wave regime of torsional shear vibrations corresponding 
to different boundary conditions at the core-crust interface and the star surface. A detailed 
analytic derivation of asymptotic spectral formulas is presented followed by a numerical 
analysis of the obtained frequency spectra. In Sec. 4, we compare frequency spectra for the 
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nodeless torsional shear oscillations computed from homogeneous and non-homogeneous 
models of the neutron star crust . The obtained results are summarized in Sec. 5. 
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2 GOVERNING EQUATIONS 

In two component model of quaking neutron star (Franco, Link & Epstein 2000), its interior 
is thought of as composed of dense core (in which self-gravity is brought to equilibrium by the 
degeneracy pressure of baryon, neutron-dominated, matter) covered by highly conducting 
metal-like crustal matter composed of nuclei (basically of iron, 26^©) dispersed in homoge- 
neous Fermi-gas of relativistic electrons. The gravitational stability of the crust is supported 
by the electron degeneracy pressure Pe related with the shear modulus as /i = n{p)pe (Blaes, 
Blandford, Madau, Koonin, 1990) with fiducial value of ratiJll /u/pe = 10~^ (Strohmayer et 
al 1991; Cutler, Ushomirsly, Link 2003). It is presumed that in the approximation of con- 
tinuous medium, the quake induced elastic deformations in crustal matter can be properly 
modeled by equation of solid mechanics for solenoidal field of material displacement Ui (e.g. 
Graff 1991; Lapwood & Usami 1981; Aki & Richards 2003) 

pu- = Vfc (Tifc (Tik = 2pUik (1) 

Uik = ^[VjMfc + VkUi] Ukk = VkUk = (2) 

expressing the second law of Newtonian elastodynamics (McDermott, Van Horn & Hansen 
1988). The restoring force is provided by elastic shear stresses aik related to shear strains 
Uik by Hooke law. Consider homogeneous model of crustal matter with constant density p 
and shear modulus /i (e.g. Epstein 1988; Bildsten & Cutler 1995; Franco, Link & Epstein 
2000). The shear character of torsional oscillations implies that they are not accompanied 
by fluctuations in the density: 6p = —pukk = 0. On substituting of atk in the equation of 
elastodynamics one has 



^ In this connection, it may be appropriate to note that the Local Density Approximation of microscopic theory of metals 
developed over the past three decades leads to the conclusion that transport coefficients of solid-mechanical elasticity, like bulk 
modulus and shear modulus /i are proportional to the pressure of degenerate Fermi-gas of conducting electrons pe squeezed 
between ions (Maruzzi, Janak & Williams 1978). The attitude that crustal matter possesses metal-like properties lends support 
to the parametrization for shear modulus as = fc(p)pe- 




(3) 



u(r, t) = a(r) a(t) 



a(t) = ao exp{iut) 



(4) 
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substitution (jl]) in ([3]) leads to the vector Helmholtz equation for the standing shear wave 

V2u + A;2u = V-u = fc^ = — ■ (5) 

In what follows all calculations are carried out in the spherical coordinates with fixed polar 
axis z. In this frame of reference the general solution of (I5l) describing axisymmetric differ- 
entially rotational oscillations of matter in the star is given by the toroidal vector field (e.g. 
Graff 1991; Lapwood & Usami 1981; Aki & Richards 2003) 

u(r, t) = V X rf/(r, t) = Vf/(r, t) x r f/(r, t) = fe{kr) P^(cos 6) exp{iujt) (6) 

f,{kr) = [Aeje{kr) + Bene{kr)] (7) 

Ur = 0, ug = 0, = fe{kr) P/(C) exp{iujt) (8) 

P^{0 = il-CY'^ C = cos^. (9) 

By je{kr) and ni{kr) are denoted the spherical Bessel and Neumann functions, respectively, 
and by P£{cos6) the Legendre polynomial of multipole degree i (Abramowitz & Stegun 
1964). Function fi obey the following recurrence relations 

f^f,_,-i±lf, f,_,,f,,,^^l±lf, (10) 

az z z 

^/£-i-(^+l)/m = (2^+l)^ z = kR (11) 

which hold for both je{kr) and n£{kr). In the long wavelength limit these functions are 
approximated by 

M^) - (^j^ Ji^iiz) = ^M^) (12) 

(2f-l)!! 2i+l 
M^) ^ —I — ni+i{z) = ne{z). (13) 

These asymptotic formulae provide a basis for obtaining analytic estimates of the frequency 
spectra of nodeless torsion oscillations. 

As a first step, we consider global torsional oscillations in the entire volume of homoge- 
neous neutron star model. For our present purpose this problem is interesting in that the 
spectral formula for the frequency of global torsional oscillations is used as the reference 
equation in testing of spectral equations for the frequency of torsional oscillations trapped 
in the crust. These latter are derived in the form showing that fiducial spectral formula for 
the global torsional oscillations is recovered when the core radius tends to zero. 

In the case of global torsional oscillations in entire volume the singular in origin solution 
of the Helmholtz equation for the displacement field must be excluded by putting Bi = 0. 
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Figure 1. The artist view of nodeless toroidal field of material displacements in the neutron star undergoing global torsional 
oscillations in quadrupole {£ = 2) and octupole {£ = 3) overtones. In quadrupole overtone (left), the field of material dis- 
placements in north and south hemispheres of neutron star undergoes out-of-phase oscillations, as pictured by arrows. In the 
octupole overtone (right), the material displacements in north and south are in one and the same phase, whereas in equatorial 
part the direction of displacements is opposite. 



Then, one has 



Ur = 0, U0 = 0, = Aiji{kr) P/(C) exp{iujt). 



(14) 



The standard boundary condition of stress free surface (7ik\j.=R = 0, where are compo- 
nents of the unit vector normal to the star surface leads to 
du, 



nr<^r(j>\r=R — 



dr 



r 







r=R 



(i-l) 



z = kR. 



(15) 



To integrate this transcendent equation one must have recourse to numerical methods. In 
what follows we focus, however, on the long wavelength limit of the last equation, when 
z << 1. This is the regime of the nodeless torsional oscillations with the field of displacements 
Ui having no nodes in the interval < r < R; in this regime of torsional oscillations the 
components of Ui are given by [ur = 0, Ug = 0, = Agr^ Pl{() exp^iut)]. From this 
the term nodeless torsional vibrations is derived. The character of shear distortions in this 
quadrupole and octupole overtones of nodeless torsional oscillations is pictured pictured in 
Fig.l (see, also, Bastrukov et al 2002). 

In the long wavelength limit, z << 1, the dispersion equation (fTSi) is reduced, with help 
of approximate formulae ([T2|) and f|T3|) , to simple algebraic relation 



z'^ - {2i + 3){i-l) = z^ 



(16) 



It follows that the frequency spectrum of global long wavelength torsional modes in a ho- 
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> 




Figure 2. Fractional frequencies v/vq, v' /vq and v" /vq of global nodeless torsional oscillations as functions of multipole 
degrees 1 < I <2Q. 



mogeneous solid star model is given by 



[{2i + 3){i-l)] uiote) 



- ^ 
P 

Un = . 

271 



(17) 
(18) 



From now on we use standard nomenclature for the torsional vibration mode, as o^^, where 
zero marks the nodeless regime of torsional oscillations of multipole degree £. 

The following equivalent representations of frequency spectrum fll8l) may be useful. First, 
given in the form 

1 



2{e + 2){e- 



2{i + 2) 



is interesting in that at large values of multipole degree, i » 1, can be replaced by 
^^i^ = [2(£ + 2)(£-l)]V2 

Second, given in the form 



i{£ + l) 



2{i + 2) 



is represented in the limit of large i as 

1^0 1^0 



(19) 



(20) 



(21) 



(22) 



The comparison of zz/z/q, v" /vq and v" jv^ is shown in Fig.2; it is implied that z/(ot^) = v. 

It is worth noting that the spectral formula (ITTll can be derived from different mathemat- 
ical footing, namely, by use of the Rayleigh's energy variational method which is particularly 
efficient when studying of non-radial nodeless oscillations of neutron stars (Bastrukov et al 
1999, 2002, 2007). 
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The spectral formulae like above derived are central to theoretical modal analysis of 
variability in electromagnetic emission of pulsars and magnetars. The works of Van Horn 
(1980) and Cioffi and Hansen (1980) were among the firsts suggested a possible association 
of microspikes of 10-50 milliseconds duration clearly discernable in the windows of main 
pulse train of radio pulsars with torsional oscillations of neutron stars (McDermott, Van 
Horn, Hansen 1988; Strohmayer 1991; Bastrukov et al 1999, 2007). This suggestion provides 
a guideline in the above mentioned current studies of quasiperiodic oscillations detected 
on the lightcurve tail of SGR 1806-20 and SGR 1900+14. The above obtained spectral 
formulae for the frequencies of nodeless torsional oscillations restored by force of shear elastic 
deformations have many features in common with those derived in (Samuelsson, Andersson 
2007) from general relativistic treatment of torsional elasticity. The detailed identification 
of QPOs detected during the flare of the above magnetars with overtones of elastic torsional 
oscillations is discussed in (Watts, Strohmayer 2007). 

3 TORSIONAL ELASTIC VIBRATIONS IN THE HOMOGENEOUS 
MODEL OF THE NEUTRON STAR CRUST 

In the reminder of this paper we focus on torsional oscillations trapped in the peripheral 
spherical layer of the neutron star of finite depth implying that seismically active zone 
depends upon the energy which is released in the starquake. The prime purpose of our 
analysis is to elucidate the effect of boundary conditions (reflecting the behavior of material 
displacements on the edges of seismogenic layer) on the form of frequency spectrum of 
elastic torsional oscillations. In so doing we adopt boundary conditions which are currently 
utilized in the works studying quake- induce oscillations in the neutron star crust. Namely 
the condition of stress-free-surface for both core-crust boundary and surface of the star and 
non-slip condition on the the core-crust interface (e.g. McDermott, van Horn, Hansen 1988; 
Strohmayer 1991; Bildsten & Ushomirsky 2000). 

3.1 No-slip boundary condition on the core-crust interface and no-stress on 
the star surface 

Torsional oscillations trapped in the crust are described by the general solution for Ui given 
by equations fl6|)- f|TT]) . To eliminate arbitrary constants Ai and Bi, two boundary conditions, 
one on the core-crust interface, at r = Rc, and second on the star surface, at r = R, must 



8 S. I. Bastrukov, H.-K. Chang, J. Takata, G.-T. Chen and I. V. Molodtsova 

be used. On the star surface we impose the standard boundary condition of the absence of 
stresses normal to surface Uk o'ik\r=R = which is reduced to 



/i 



dr 



r 







J r=R 



and on the core-crust interface the no-shp condition 



(23) 



(24) 



Physically, this condition means that the amplitude of differentially rotational oscillations 
triggered by starquake in the neutron star crust is gradually depreciated from the star surface 
to the core-crust interface which is the internal boundary of seismogenic layer. On inserting 
([121) in (El and (El we obtain 



fe+i{z) [i- 







Xz 



z = kR 
^ A < 1. 



(25) 
(26) 



Note, the notation Zc = Xz means that the radius of the core Rc can be represented as 
Rc = XR with A from the interval ^ A < 1. It is worth emphasizing that A is strongly less 
than unit. In terms of Bessel and Neumann functions these latter boundary conditions read 



Je+i[z) 



(£-1) 



ni+i{z) 



nt{z) 







Aat{Xz) + Btnt{Xz) = 0. 



(27) 
(28) 



Casting these equations as homogeneous matrix equation whose Wronskian must be equal 
zero we arrive at the dispersion equation of the form 



W{z) 



ne{Xz) 



nt+i{z) 



uAz) 



0. 



(29) 



Computations of roots of this transcendent equation in which spherical Bessel and Neumann 
functions are defined in different points is non-trivial numerical problem (e.g. Pexton & 
Stiger 1977). However, our prime purpose here is to evaluate the frequency spectrum for long 
wavelength differentially rotational oscillations, when z = kR << 1. Taking into account 
that in the long wavelength limit 

-1 (2^-1)!! 



Je 



{z)ne{Xz) P A-(^+^) UXz)n,{z) -^(3X' (3 = -z' 



(2£ + l)!! 



one finds that exact dispersion equation (1291) is reduced to 



= (2i + 3)[(£ -!) + {£ + 2)A2^+i] z^ = k^R^ = ^R\ 



(30) 



(31) 
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From this it follows 



= (2£ + 3)(£-l) 



1 + 



- 1 



1^1 



— = 1 — n. n = . 

R R 

In the limit of zero-size radius of the core, A = {Rc/R) 



UJl UJQ 

— z/0 = — A 

ZTT ZTT 



(32) 
(33) 



0, corresponding to torsional 
oscillations in the entire volume of a solid star, we regain spectral equation (ITS]) for the 
frequency of global torsional mode. For our further purpose we note that equation (!32|) can 
be represented in the following equivalent form 



^0 



2 1- 



1 



1 + ^±l;,2^+i 



1/2 



(34) 



2(£ + 2); V £-1 

which is discussed in the next subsection. Henceforth the suffice in expressions for v and uj 
marks the number of eigenfrequency problem under consideration. The considered in this 
subsection is marked by suffice 1, and two another conceivable boundary conditions are 
regarded in the reminder of this work. 



3.2 Stress free boundary conditions on both core-crust interface and the 
neutron star surface 

Now we adopt the boundary conditions of the free from stresses surfaces for both the surface 
of the star and the core-crust boundary aik\r=R,R^ = 0: 

duA, u„ 



dr 



nrar^\r=R,Rc — 

J r=R,Rc 

which in terms of fe{z) are given by 







dfejz) fejz) 
dz z 







dZr 



^ A < 1. 



(35) 



(36) 



- Zc=XkR 



For a non- vanishing solution of these equations to exist, the following dispersion equation 



Je+i{z) (« - 1) 



je+i{Xz) - 



\z 



n,+i{\z) - ^^{i - I) 
Xz 



, , nAz) ^ , 
n,+^{z)-^^{^-l) 







(37) 



must hold. Following the line of argument of foregoing subsection, consider the limit of long 
wavelengths, z = kR << 1. This yields 



z/, 



4 = (2£ + 3)(£- 1) 



1 + A^^+^ 
1 - A2^+3 



^ A < 1. 



(38) 



In the limit A = {Rc/R) —>■ 0, the last equation is reduced to spectral formula f|T8|) for the 
frequency of global torsional oscillations. 
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k: 1 1 1 1 

15 10 15 20 



/ 

Figure 3. Fractional frequencies of nodeless torsional shear oscillations as functions of multipole degrees 1 < £ < 20 plotted at 
h = AR/R = 0.05,0.1 and 0.3. u/uq for global nodeless torsional mode in entire volume of homogeneous neutron star model; 
fi I vq and 1^2 / I'D for the frequency of torsional modes trapped in the peripheral layer of homogenous model of crust computed 
with boundary conditions of first and second examples. 



Fig.3 illustrates the general trends of fractional frequencies - zz/z/q, z/i/z/q and z/2/z/o as 
functions of multipole degrees 1 < £ < 20 plotted at the values of fractional depth given 
by = Ai?/i? = 0.05, 0.1 and 0.3. The difference caused by boundary conditions is notable 
in the low-^ domain, when the depth of seismogenic layer A/2 = 0.5 and 1.0 km (in the 
neutron star of fiducial radius i? = 10 km). This difference is practically vanishes when the 
thickness of layer is about Ai? = 3.0 km and larger. Generally, the thicker this depth the 
less difference. Also we note, for our further purpose, that equation ( l38l) can be represented 
as 



^ = [(£ + 2)(^-l)]V^p,-i 



2 1 



2/ l-A2^+3 



(39) 



At this point it seems worthy of noting that analogous, from mathematical side, problem 
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of torsional vibration mode trapped in the peripheral solid layer has been considered long 
ago by Pekeris (1965) in the geoseismic context. The obtained in this latter work spectral 
formula with boundary conditions identical to considered in this subsection reads 

^ = [(£-l)(£ + 2)]V2p-i £^10. (40) 

In work of Pekeris (1965) the dimensionless parameter p£, has been numerically computed 
from integral equation establishing "relation between geometrical optics and terrestrial spec- 
troscopy" and tabulated for 1 ^ £ ^ 50. Remarkably, in appearance the Pekeris spectral 
formula is identical to the above presented equations (134|) and (139|) . In Fig.4 we plot t'p/z/Q 
with use of data for pi taken from (Pekeris 1965, Table 2) as a function of i in juxtaposition 
with our spectral formula (HUj) in which we have used the value of fractional depth of seis- 
mogenic layer given hj h = 0.1. Analogous conclusion holds for vi/vq. It should be noted, 
however, that equation pUj) has been obtained with use of non-uniform profile for the speed 
of transverse wave of elastic shear Ct(r) which is used as input function of the method, not 
computed. In view of this, the last figure can be considered no more than juxtaposition, 
not a comparison, of our and Pekeris spectral equations. Nonetheless, the fact that both 
approaches yields practically identical results suggests that presented in this section analy- 
sis can be utilized for assessing frequencies of torsional seismic vibration modes in the solid 
Earth-like planets too. 
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4 NODELESS TORSIONAL OSCILLATIONS IN HOMOGENEOUS AND 
NON-HOMOGENEOUS MODELS OF THE NEUTRON STAR CRUST 

From above it follows that in the long wavelength limit the toroidal field of displacements 
(!6l)- (fTT]) is reduced to the form that can be conveniently represented as follows 

u(r,t) = a(r)a(t) a{t) = aoexp^iut) (41) 

a(r)=Vx(r)xr x{r) = Mr)Pe{0 (42) 

Mr) = [Aer' + Ber-^''-% (43) 

The toroidal field (HTi) -( l43l) is the general solution to the vector Laplace equation 

V^u(r,t)=0 Vu(r,t) = 0. (44) 

which can be thought of as the long wavelength limit of the Helmholtz equation V^u + Zc^u = 
0, because in the limit of long wavelengths, A oo, the wave vector k = (27r/A) —>■ 0. The 
vector Laplace equation can be regarded as fundamental equation defining regime of nodeless 
shear vibrations q The practical usefulness of this attitude is that it allows to assess the 
difference between predictions for the frequency spectra computed on equal footing within 
homogeneous and inhomogeneous models of the neutron star crust, that is, regardless of the 
form of density and shear modulus profiles. Most efficiently it can be done by Rayleigh's 
energy method (Bastrukov et al 2007). The point of departure is the integral equation of 
energy conservation which is obtained from equation of elastodynamics ([T]) by its scalar 
multiplication with Ui and integration over the crus volume 

— J — = ^ J ^ikUik dV = -2 J jiUikiiikdV. (45) 

On inserting in fj45l) . the separable form of the displacement field f HTj) . we arrive at equa- 
tion for temporal amplitude a{t) having the form of the well-familiar equation of normal 

vibrations 

dE ^ Ma^ Ka'^ 

E = — — + — - (46) 



dt 2 2 

K 

M 



K 

a + u'^a = = 77 (47) 



M 



j p{r)aiaidV K = 2 j p{r) aik aik dV a^fc = ^[VjOfc + VfeOi]. (48) 



2 



The axisymmetric, odd parity, toroidal vector field represents one of two fundamental (mutually orthogonal and different 
in parity) solutions to I I41I I describing nodeless torsional vibrations. The second fundamental solution is given by even parity 
poloidal vector field (Bastrukov et al 2007) which describe nodeless spheroidal vibration mode, in accord with canonical Lamb's 
classification of vibrational modes in an elastically deformable solid sphere (e.g. Lapwood & Usami 1981; Aki & Richards 2003). 
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The analytic form of the inertia M and stiffness K shows that method can be utihzed for 
computing frequency u"^ = K/M of shear vibrations (both spheroidal and torsional) for wide 
class of models with non-uniform density and shear modulus profiles in the crust which are 
the input parameters of the method. Taking the integral over the solid angle in the above 
expression for inertia M we obtain 



M 



2i + l 



R 



p{r) r'^^'^'^dr + 



Ai 



R 



R 



^2 

p(r) rdr + — / p(r) r~'^'^dr 
Ai 



Rc Rc R> 

In similar fashion, for the parameter of rigidity K we get 



(49) 



K = A7TAiit-l 



R 



jj,{r) r dr + 



+ 2) 



R 



A\ (£ 



/i(r) r 



-2i~2 



dr 



(50) 



These later equations for M and K emphasizes the fact that the choice of boundary condi- 
tions does matter for the problem under consideration. 

As a representative example, consider torsional oscillations in the crust with constants 
Ae and Be eliminated from the following boundary conditions 



U^\r=Rc = U^\r=R = [4>R X R-]. 

4)ji = exp{iujt)VnPi{() Vn = 



0, 



d 



d 



(51) 
(52) 



89 sin 9 d(j) ) 

First is the no-slip condition on the core-crust interface, r = Rc-, and the form of second 
boundary condition on the star surface, at r = i?, is dictated by symmetry of general 
toroidal field of nodeless torsional oscillations. The resultant algebraic equations steaming 
from boundary conditions fl^ lead to 



A, = U, Bf = -UfRf+^ 



^21+1 _ ^21+1 ■ 

In the homogeneous crust model, presuming constant values of p and p, we get 

47Ti{i + 1) pR^ 



(53) 



M{i, A) 



(2£ + l)(2£ + 3) (l-A2^+i)^ 



X 



(2£ + 3)A 



2e+i 



+ 



{2i 



4ni{e - 1) 



2i-l 
pR^ 



2i+l (1 - A2^+i 



^2U3 
1 - 



2i + 3 
21-1 



y2(2l+l) 



+ 2) y2l+l 



(54) 
(55) 



(£-i: 

Note, in the limit A — ^ we again recover the spectral formula for the frequency uo"^ = K/M 
of the global torsional oscillations in the entire volume of the homogeneous solid star model 

m 
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/■[10 km] 



Figure 5. Density profile for realistic non-homogeneous neutron star model built on the TOV equation and EOS given in 
(Wiringa ot al 1988). Insertion depicts density distribution in the peripheral crustal region. 
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Figure 6. Profile of speed of transverse shear wave in the crust of non-homogeneous model with density pictured in Fig. 4. 



4.1 Numerical analysis 

In the above expanded the energy variational method of computing the frequency of elastic 
shear oscillations, the profiles of density p(r) and pressure p(r) linearly proportional to the 
shear modulus /i(r) of crustal matter are regarded the input parameters. The approach we 
adopt is the one used by previous authors taking these parameters from evolution models 
aimed at computing gravitationally equilibrium state of matter and the internal structure 
of neutron star. In these later models the density and pressure are computed from Tolman- 
Oppenheimer-Volkoff (TOV) equation with account for realistic EOS (e.g. Weber 1999; 
Lattimer & Prakash 2001). 

Specifically, we adopt here the model of neutron star model with fiducial mass of M = 
1.43Ms„„ presented in Wiringa, Fiks & Fabrocini (1988) and will use parameters for the 
crust matter given in Douchin & Haensel (2001). The density profiles in the star as a whole 
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and in the crust of this star model are pictured in Fig. 5. The position of the core-crust 
interface is placed at the fiducial density p = 1.5 10^^ g cm~^, in accord with arguments of 
works (Cutler, Ushomirsky, Link 2003; Pethick, Ranenhall & Lorentz 1995). The effect of 
non- homogeneous density in the crust on the speed of transverse wave of elastic shear in 
this model is illustrated in Fig.6. 

In Fig. 7 we compare the fractional frequency of torsional oscillations trapped in the 
crust of fixed depth AR = 0.6 km (that has been normalized to the value of uq = 15 Hz) 
computed in the homogeneous and inhomogeneous crust models. It is seen that prediction 
of homogenous crust model with boundary conditions of this section are quite different from 
that inferred in previous section from homogeneous model too, but for different boundary 
conditions. Also, the predictions of homogeneous models are drastically different from those 
for non- homogenous one. This difference is manifested in both the overall trends of frequency 
as a function of multipole degree and absolute values of fractional frequencies. One of sources 
of uncertainties is parametrization of shear modulus profile. In the curve computed with 
fi = k{p)pe we used parametrization of work (Cutler, Ushomirsky, Link 2003) and with 
fi = cfp from (Strohmayer et al 1991). At this point we leave the discussion of mathematical 
details of numerical calculations because this issue is not the main subject of presented 
investigation. 
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5 SUMMARY 



There are now quite solid arguments showing that the Soft Gamma-Ray Repeaters are iso- 
lated, non- accreting, seismically active magnetars (Thompson & Duncan 1995) - quaking 
neutron stars endowed with ultrastrong magnetic fields. The X-ray bursting luminosity of 
magnetar is associated with starquakes which are thought of as sudden release of magnetic 
field stresses cracking of the neutron star crust by the X-ray fiares. Evidence in favor of seis- 
mic nature of the magnetar flares is provided by striking similarities between statistics of 
SGR's bursts and earthquakes (Cheng, Epstein, Guyer & Young 1996). The common behef 
is that the dynamics of quake induced internal elastic and magnetic field stresses in the crust 
can be properly understood within the framework of the two-component core-crust model of 
quaking neutron stars (Franco, Link & Epstein 2000), provided that the core-crust couphng 
is dominated by an ultrastrong magnetic field. A conceivable and comprehensive, from a 
physical point of view, explanation of the magnetic core-crust coupling provides a model of 
paramagnetic neutron star (Bastrukov et al 2002; Bastrukov et al 2003). In this model the 
neutron star core is regarded as a spherical bar magnet (composed of baryon matter domi- 
nated by poorly conducting neutron component) permanently magnetized to saturation due 
to Pauli's mechanism of field-induced paramagnetic spin-polarization of neutron magnetic 
moments. The microstructure and extremely high conductivity of crustal matter indicates to 
its metal-like electrodynamical properties. This analogy suggests that the magnetic coupling 
of the core (permanent magnet) with the crust (metal) is similar to the well-known magnetic 
adhesion between a metal specimen and bar magnet. The fast process of postquake recovery 
of the magnetar is of course dominated by force of gravitational pull. The prime effect of 
magnetic field frozen in the core on this process is that the the perturbed by starquake 
a highly conducting crustal matter (as well as the less dense plasma of magnetar corona 
(Beloborodov & Thompson 2007) expelled from the star surface by SGR's flares) sets in ax- 
isymmetric torsional oscillations about axis of ultrastrong magnetic field frozen in the star. 
It is these oscillations about magnetic axis are observed, as is argued in (Israel 2007; Watts, 
Strohmayer 2007) as quasiperiodic oscillations of x-ray luminosity with high signal-to-noise 
ratio. In (Watts & Strohmayer 2007) arguments are given that the data favor the idea of 
a seismic origin of the detected QPOs, that is, as caused by quake induced torsional shear 
oscillations of crustal matter of magnetar. Adhering to this attitude, we have investigated 
several scenario of torsional shear vibrations restored by bulk force of shear elasticity (the 
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problem from which the very notion of torsion shear vibrations of an elastic sphere came 
into existence). Our prime purpose was to elucidate the distinctions between spectral formu- 
lae for the frequency of nodeless torsional shear oscillations caused by different behavior of 
quake induced crustal matter on the boundaries of the crust. In pursuing this goal and using 
canonical methods of mathematical physics we have considered several examples of exact 
solutions of the eigenfrequency problem demonstrating crucial effect of boundary conditions 
on the frequency spectrum of torsional mode of elastic oscillations trapped in the crust. The 
asymptotic spectral formulae for the frequency of nodeless torsion oscillations have been 
presented so that they can be conveniently applied to a wide class of celestial objects. The 
practical usefulness of the exact solutions for the toroidal field of material displacements con- 
sidered here is that they can be utilized in the study of torsional shear vibrations restored 
by forces of intrinsic stresses of different physical nature, like Newtonian gravitation field 
stresses (e.g. Shu 1992) and Maxwellian magnetic field stresses (e.g. Chandrasckhar 1961; 
Mestel 1999), not only Hookean elastic one. With all that we conclude, while the obtained 
spectral formulae properly match the observable QPOs frequency, the information obtained 
from only these models is less than necessary to draw definite physical statements which of 
above forces plays dominant role. The point of particular interest, as is generally believed, 
are torsion oscillations driven by the Lorentz force that can be represented as divergence of 
fluctuating magnetic field stresses (Franco, Link, Epstein 2000) 

5Tik = -;—[Bi 5Bk + Bk 5Bi — Bj SBjSik] 

4:71 

= V X [u X B] V • u = 0. 

The canonical form of magneto-solid-mechanical equations governing Alfvenic, compression 
free, oscillations in a perfectly conducting elastically deformable sohd, regarded as material 
continuum, reads 

The integrand of equation of energy conservation exhibits the fact that the work done by 
magnetic field stresses in the volume of a quaking neutron star is accompanied too by 
shear deformations and, thus, shows the torsional shear oscillations can also be sustained by 
fluctuations of the magnetic field stresses. We shall turn to this problem in a forthcoming 
paper. 

The authors are indebted to Dr. Judith Bunder (UNSW, Sydney) for critical reading of 
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